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Abstract 

A manifestly covariant expression for the current matrix elements of 
three quark bound systems is derived in the framework of the Point 
Form Relativistic Hamiltonian Dynamics. The relativistic impulse 
approximation is assumed in the model. A critical comparison is made 
with other expressions usually given in the literature. 

PACS number(s): ll.30.Cp, 24.10.Jv, 13.40.Gp 
1. Introduction 

Aim of the present work is to show that, in the context of the Point Form 
Relativistic Hamiltonian Dynamics (PF RHD), it is possible to write, for 
hadronic bound systems, manifestly covariant matrix elements of the current 
operators. 

For clarity, we define as manifest covariance the property of an equation of 
being written in terms of quantities that (i) transform in a well-known way 
under Lorentz transformations and (ii) are not related to a specific reference 
frame. 

We use a relativistic impulse approximation [1,2] (RIA), that generalizes at 
a relativistic level the widely known model used for the study of the electro- 
magnetic interactions of nuclear systems. 

We shall examine in more detail the case of the four-vector electromagnetic 
current, but the method has been developed to study any kind of interaction. 
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in particular the axial-vector case [3,4] that is relevant for the weak structure 
of the hadrons. 

In this article we shall specifically refer to the nucleon as a system composed 
hy N — 3 (identical), s — 1/2, constituent quarks. The same procedure can 
be immediately generalized to the study of the N ^ 3 composite systems. 
The cases of nonidentical constituents and s ^ 1/2 will be studied, within 
the same theoretical framework, in subsequent works. 

We also show that our model, that is the result of an independent investiga- 
tion [5], is completely equivalent to the standard formalism (SF) of PF RHD, 
developed in refs.[l-3] and implemented with great clarity and precision in 
ref.[6]. In this concern we anticipate that the main differences between our 
model and SF are the following: 

(a) we use, in our work, the formalism of the Dirac equation boosting covari- 
antly the Dirac spinors, while in SF Wigner rotation matrices (not manifestly 

covariant) are employed; 

(b) the spatial part of the bound system four-velocity and the (independent) 
three-momenta of two quarks are used in our work as spatial variables for the 
representation (or projection) states, while the velocity states are considered 
in SF; the use of these states, that are related to the rest frame of the bound 
system does not fulfill the requirement (ii) given above to have manifestly 
covariant equations; 

(c) furthermore, we perform manifestly covariant integrations over that spa- 
tial variables to calculate the matrix elements of the current operators. 

We highlight that the electroweak matrix elements calculated by means of 
PF RHD allow to reproduce with good accuracy the experimental nucleon 
form factors [2-7]. 

The model studied in the present work allows for two further developments 

(to be studied in different works): the definition of a dynamically conserved 
electromagnetic current by means of a suitable minimal coupling procedure 
and the introduction of (virtual) negative energy states in the current matrix 
elements. 

The paper is organized as follows. 

In sect. 2 we revise the construction of the Poincare algebra generators, in- 
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troducing, at the same time, the operators that are used to describe the 

dynamics of the bound system. 

In sect. 3, by introducing the projection states of our model, we define the 
wave functions for the bound system in the framework of PF RHD, also 
discussing their boost properties. 

In sect. 4 we explicitly construct the Dirac wave functions of our formalism, 
showing the equivalence of their boost properties with those of the SF. 
Finally, in sect. 5, the matrix elements of the current operators are studied by 
means of the RIA. The main result of this work is our covariant expression 
given in eq.(5.3c). An accurate comparison with SF is performed transform- 
ing eq.(5.3c) into the standard form of eq.(5.20). 

2. The Poincare Algebra 

In the present work, considering particles of mass m, always on-shell, we 
transform the four- momentum p'^ — (e(p), p) by means of a canonical boost 
written in the following standard form 

e6(p; v) = e(p6(p; v)) = v°e{p) + vp (2.1a) 
Pb{p; v) = p + v( vp^^ + e(p)) (2.16) 
The two previous equations are usually resumed in the form 

In eqs.(2.1a,b) we have introduced the time component of the four-momentum 
of the particle, i.e. the energy, as e(p) = [p^ + m^]^/^ and = ([v^ + 1]^''^, v), 
that is the four- velocity boost parameter. We recall that the physical veloc- 
ity of the initial frame measured from the boosted one is u — v/v^. 
The independent transformation equation, that is used to define the boost 
in the Hilbert space, is eq.(2.1b), while eq.(2.1a) can be obtained from that 
one by calculating the on-shell energy of the particle with the boosted mo- 
mentum. 

As anticipated, we follow the scheeme of the PF RHD when defining the 
generators of the Poincare algebra [1]. In more detail, for a system of three 
quarks, the total angular momentum J and the total boost K, being free of 
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the interaction, are written as the sum of the single particle generators, in 
the form: 



J = ^ (r, X p, + s,) (2.2a) 



i=l 



(2.26) 



2 e(pi)+m_ 

where pj, r^, Sj = l^j, m and e(pj) respectively represent the three-momentum, 
the conjugated (position) variable, the spin, the mass and the energy of the 
i-th quark. 

For completeness we also give the expression of the finite boost operator, 
that in the PF RHD is not modified by the interaction: 

5(v) =exp(iK-U) ~ l + i5u-K (2.2c) 

with 

U = ^ tanh-^(^) = tanh-^(|u|) (2.2(i) 



On the other hand, the total four-momentum operator of the system, that is 
pt^ = (po = H, P) depends on the interaction among the constituent quarks. 
We shall define the operator P^^ in eq.(2.19). To this aim, we have, previously, 
to introduce: 

(a) the quantum-mechanical operator V^, that represents the four- velocity 
of the bound system measured from a generic reference frame (GF); 

(b) the other dynamical variables of the quantum mechanical model. 

We first consider point (a), that is the construction of the quantum me- 
chanical operator V^. In order to help the reader to understand the physical 
meaning of the following procedure, we note that the four-momentum P'^ of 
a system, as a classical quantity, can be written in terms of as 

P^' = MV^ (2.3) 

where the physical mass M of the (bound) system has been introduced. The 
corresponding quantum-mechanical expression will be given in eq.(2.19). 
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To derive this expression, we have to write as a function of the momenta 
of the constituents. As first step, we introduce the rest frame (RF) four- 
momentum of the i-th quark 

pr^{e{p*),P*) (2.4) 

Here and in the foUowing, the asterisk denotes the quantities observed in the 
RF. The sum of the four- momenta of the three constituent quarks, is, by 
definition of the RF 

Epr = (Ee(pn = M,,o) (2.5) 

i=l i=l 

where we have also introduced Mf that represents the free mass operator 
of the system. By applying the Lorentz transformation of eqs.(2.1a,b) (as 
a function of the parameter V) to the p*'^ and also using eq.(2.5), one can 
write the sum of the four-momenta of the particles in a GF as 

= V^Mf (2.6) 

i=l 



with 



Pi = (e(Pi),Pi) 



We highlight that Mf, as defined in eq.(2.5), is a nonvanishing and Lorentz 
invariant quantity. The nonvanishing character of Mf allows to solve the 
previous equation with respect to V^. Lorentz invariance allows to write Mf 
in terms of the observed in a GF. In this way one can express as a 
function the p^, or, more precisely, of the three-momenta pj. 
(Note that, the pj, with i = 1,2,3 represent, in the first step of the con- 
struction, the spatial dynamical variables of the relativistic model. The final 
choice of the spatial variables will be given in the following.) 
In more detail, Mf is expressed as a function of the momenta in a GF in the 
form 



Mf = M/(pi,p2,P3) 



3 

ij=l 



-I 1/2 

(2.7) 



that will be taken as the definition of the operator Mf. In consequence,we 
can also write 
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^'^(Pl,P2,P3) = [M^(pi,P2,P3)]-'EPi 



(2.8a) 



1=1 



and , obviously 
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(2.86) 
(2.8c) 



Let us note that the observable four- vector V^, as given in eq.(2.8a), trans- 
forms in the same way as a standard four-momentum, that is replacing p 
with V, e(p) with y°(V) in eq.(2.1b). In this way we introduce 



This result, that is also consistent with eq.(2.3), can be easily derived by 
transforming, with the help of eqs.(2.1a,b), the that appear in eq.(2.8a). 

As for point (b), we can now introduce the final choice for the complete 
set of commuting operators that will be used for the quantum mechanical 
description of the system. To this aim we note that, due to its definition in 
eq.(2.8a), the operator commutes with the momenta of all the particles. 
In consequence, it is possible to choose the following operators: 

(i) as spatial variables, the three-momenta of 2 quarks, say P2,P3, and the 
spatial components of the four-velocity V; those variables replace the first 
step choice of Pi, P2, Ps; 

(ii) the spin operators of the three quarks; the eigenvalues of their projections 
on the z axis will be denoted as (7i, (72, (73. 

For further developments, it is necessary to express Pi and Mf as functions 
of P2, P3 and V. 

First, we recall that the rest frame quark energies are invariant quantities 
[8] , that can be written as 



n° = n°(V;v)=y'^(V,(V;v)) 



(2.9a) 



V6 = V5(V;v) 



(2.96) 



(2.10) 



Second, we write eq.(2.6) in the form 



+ = -P'l + ■ [e(pt) + 6(p;) + 6(p*)] 



(2.11) 
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Then, squaring both sides, with the help of eq.(2.10), one obtains the RF 
energy of the quark #1 as a function of pi, p2 and V: 



e(pD = et(p2,P3, V) = 

- (P'2+P'3)(P2+P3)9,. + [(p'2+p's)V,]f^' (2.12) 

where P2, and are functions of p2, Pa, and V, respectively. 
Finally, we find 

M/(p2, P3, V) = (p^ + p^,)V, + el{p2, P3, V) (2.13) 
and, by means of eq.(6) 

(P2, P3, V) = -(p^ + p^) + ■ Mj(p2, P3, V) (2.14) 

By definition, Mf is a Lorentz invariant operator, that is 

[K,Mf] = (2.15) 

We now introduce the interaction among the quarks by means of the total 
mass operator M that, according to the Bakamjian- Thomas construction 
[1,9], is defined as 

M = Mf + W (2.16) 
where W reprents a Lorentz invariant interaction operator, that means 

[K, W]^0 (2.17a) 

and, in consequence, 

[K, M] = (2.176) 

In this work, we do not enter into the details of the definition of W. We 
only point out that rotationally scalar operators, defined in the RF (as the 
phenomcnological potentials generally used for the relativized constituent 
quark models, in particular the hypercentral potentials [10]), are formally 
Lorentz invariant and can be also written in an explicit invariant form by 
means of the dynamical variables of the model. 

Note that, if the interaction operator W represents a quasi-potential derived 
from an underlying field theory, its expression is, in general, highly momen- 
tum dependent. 
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In any case, being defined in the RF, the interaction operator W has nonva- 
nishing matrix elements only between states with the same V^, that means 

[V^',W] = (2.18) 

We can now introduce the generators of the time and space translation of 
the system, that is the four-momentum operator, as 

pi^ ^ M ^ (Mf + W) ■ V (2.19) 

that is the same expression of eq.(2.3), but considered as a definition of a 
quantum mechanical operator. 

Standard calculations [1] show that the total generators defined in eqs.(2.2a,b) 
and (2.19) fulfill the Poincare group commutation rules [11,12]. 

3. The wave functions of the model 

We now turn to introduce the representation states that will be used to write 
down explicitly the wave functions of the model. Following the definitions of 
the dynamical variables given in the previous section, one has 

IV'r >= |P2, Pa, V; C7i, C72, (73 > (3.1) 

with the normalization 

< '^^rWr > = < P2, P3, V; CTi, CTs, CTglPa, P3, V; a[, CT2, CT3 > = 

S{P2 - P2)S{Ps - P3)S{y' - V)(5o-i'a/(T2'(T2(^(T3'o-3 (3-2) 

The choice of eq.(3.1), as it will be shown in sect.4, helps to introduce in a 
very clear way the relativistic impulse approximation for the current matrix 
elements. On the other hand, in SF a difi'erent type of representation states, 
currently denoted as velocity states, is generally used to study the relativistic 
bound state wave functions. In the velocity states the spatial variables are 
represented by V and by the three (not indcpcdcndcnt) rest frame momenta 
P15P25P3 better by the two (independent) Jacobi momenta Pp,PA- As 
shown in ref.[2], the Lorentz transformation of these states is given by the 
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standard boost of V and by a Wigner rotation of the rest frame momenta or 
of the Jacobi momenta. If also the spin projections are referred to the RF, 
the same Wigner rotation acts on the spin variables. 

By applying the boost operator of eq.(2.2c) to the representation states, one 
obtains: 

-B(v)|p2,P3,V; (71,(72,(73 > = 
-R(Pi;v)-R(p2;v)i?(p3;v)|p26,P36,V6; (71,(72,(73 > G(p2, P3, V; v) (3.3a) 
with 

G(P2,P3, V,V) = [ ,(p^) ,(p3) VO(V) ] 

The previous equations show that the action of the boost operator on the 
representation states can be divided into a spatial (a) and a spin (b) part. 

(a) The spatial part, denoted in the following as -B(v), produces an eigen- 
state of the boosted momenta p26, P36, that are taken as functions of the 
corresponding unboosted variables by means of eq.(2.1b); the numerical fac- 
tor G(p2,P3, V) is due to the nonhnearity, with respect to the momenta, of 
the boost generator of eq.(2.2b) and provides for the correct normalization of 
the state, being -B(v) a unitary operator. To simplify further developments 
we introduce the following spatial matrix element 

< P2, P3, V|B(v) IP2', ps', V >= G'(P2, P3, V; v) 
^^(P2 - P26(P2'; V))5(P3 - P36(P3'; v))5(V - V,(V'; v)) (3.4) 
and recall the following property of the delta functions 

KVi - Pi6(p/; v)) = 5(p/ - Pi6(pi; -v))-^%^ (3.5) 

e&(Pi ; v) 

with i = 2, 3. Note that Pi6(Pi; — v) represents the inverse Lorentz transfor- 
mation on pj that is obtained using in eq.(l.b) the boost parameter —v. 

(b) In eq.(3.3a) the spin part is given by the product of the i?(pj; v) that 
represent the Wigner spin rotation operators (due to the second term in 
the generator of eq.(2.2b)) that depend on the numerical values of the pj. 
The (not independent) momentum pi is obtained by means of eq.(2.14). By 
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considering the Pauli spinor representation for the spin states, for further 
developments we introduce the following matrix elements 

W^'RiPh v)Wa, = Ra^'a, {PU v) (3.6) 

In SF the matrix elements of the spin rotation operators have been denoted 
as 

i?..v,(Pi; v) = Dl!;^jRw{Pi,B{v))] (3.7) 

Such notation is used to represent the spin 1/2 rotation matrices considered 
as functions of the Wigner rotation related to the momentum pi and to the 
boost -B(v). 

The wave function of our model is determined in the RF, as a function of the 
Jacobi momenta Pp, Pa solving the mass eigenvalue equation for the mass 
operator introduced in cq.(2.16). This sohition is a velocity state solution 
with V = 0, as indicated in the next equation by a Dirac S function. It is 
written as 

^rf{p„ Pa, V) = ^l^'^'ip,, px)5{Y) (3.8a) 

with 

V''^^(Pp,Pa) = Yl V'aia2(T3(Pp'PA)«'ai'ti'<T2li^<T3 (3.86) 

lTl(T2Cr3 

In the previous expression J, E respectively represent the total angular mo- 
mentum (absohitc value) and its projection on the z axis. This state is con- 
structed by coupling the angular momenta with Clebsch-Gordan coefficients, 
for example according to the standard scheeme [13] 

with 

Note that, in eq.(3.8b) the dependence on the quark Pauh spinors has been 
highlighted in order to make a comparison with SF. 

For the following developments, it is convenient to introduce, as spatial vari- 
ables, instead of the Jacobi momenta, the RF momenta P2, P3. The former 
and the latter momenta are connected by a standard linear relation. We have 
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^'''{P2, P3) = j'/'V''''^(Pp(p;, Pi), Px{P2, P3)) (3.9) 

where j'^^^ represents the (numerical) constant factor that is used to keep the 
normahzation to unity for the wave function when using the new variables 
P2, P3. The wave function of the previous equation can be decomposed with 
respect to the Pauli spinors in the same way as the wave function given in 
eq.(3.8b). As before, 

^rf{p;, P3, V) = i^'^ipl, p;)5(V) (3.10) 

We can now determine the wave function of the system in a GF boosting 
the RF wave function given in the previous equation. We use the boost 
parameter vg, that, as usual, represents the spatial part of the four-velocity 
of the bound system observed from the GF. One has 

V'g(P2,P3,V) =< P2,P3, V|S(vg)|V'RF > = 
J d^p^'d^Ps'd^y < P2,P3, V|B(vg)|P2',P3', V > 

<P2,P3,y'\^RF> (3.11) 

By using the explicit expression of the RF wave function of eq.(3.10), the 

property of the spatial part of the boost and of the S functions, respectively 
given in eqs.(3.4) and (3.5) and, finally, the spin rotation operators of eq.(3.6), 
one obtains 

V'g(P2, P3, V) = R{pIg', vg)^(P2g; vg)^(P3g; vg) 

<P2,P3,V\B{vgMrf> (3.12a) 
where we have introduced the spatial part of the boosted wave function 

"e(P2G)e(P3G)^'^' 



< P2,P3, V|5(vG)|V'i?F > = 



:i+v^)-v4 



e(P2)e(P3) 

^^'^(p;G,Py'^(V-VG) (3.126) 
also, by means eq.(2.1b), we have used 

P*G = PbiPh -Vg) (3.12c) 
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that represent the rest frame three-momenta considered as functions of the 
three-momenta of the GF, transformed by means of the parameter vq. 

4. The Dirac equation formalism 

In order to construct operators that manifestly transform as Lorentz tensors, 
it is very useful to make use of the Dirac equation formalism. 
First, we define the RF Dirac wave function in the form 

Kf(P2, pI V) = u(pI)u(p*)u{pI)^Ijrf{p*, pI V) (4.1a) 
with the Dirac spinors, 

^e{pi)+m 



For brevity, we denote these quantities, here and in the following, as (positive 
energy) Dirac spinors, taking into account that they represent 4x2 matrices 
acting onto the Pauli spinors contained in ■^iii? ( P3, V). They are 
covariantly normalized as u{pi)u{pi) = 1. 

We recall that the Dirac spinors are boosted by means of the nonunitary 
Dirac boost operator 

5f (v) = [Sf (v)]+ = [i(.° + 1)]V^ + [1(.° - i)]V2 (ZO^ (4.2) 

where we have introduced the Dirac the gamma matrices '-fi'^ = (7i°,7i) for 
the i-th particle; also, is the time component of the four- velocity boost pa- 
rameter. Standard calculations show the following very important property 
of the Dirac boost when applied to the Dirac spinors 

BfiYHp,) = 7i(p5(p.; v))i?(p,; v) (4.3) 

It shows that the Dirac boost produces a Dirac spinor of the boosted mo- 
mentum applied to the spin rotation operator, given in eq.(3.7), that acts 
onto the Pauli spinor. 

Introducing 

S^(vg) = 5f (vg) ® B^{yg) ® B^{yg) (4.4) 
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we now construct the GF Dirac wave function for the three quark system by 
means of the following boost 



^g(P2,P3, V) = B^iYcHplaHplaHpla) 

< P2,P3, V|5(vg)|V'rf >= (4.5a) 

= u{Pl)u{P2)u{p3)lpG{P2, P3, V) (4.56) 

where eqs.(3.12a,b) and (4.3) have been taken into account. Also, equiva- 
lently, making explicit use of eq.(3.12b), one can write 

1/2 

<^g(P2, Pa; vg)(1 + v^)-V45(V - vg) 

(4.5c) 
with 

<^g(P2, Pa; vg) = S^(vg)m(pIg)w(P2g)«(P3g)^^''(P2G, P3g) (4-5c^) 

The expression V'g (P2? PS) ^) of eq.(4.5c) is the boosted Dirac wave function 
of the model. Also, (fioiP^, Pa; vg) of eq.(4.5d) can be defined as the boosted 
intrisic Dirac wave fuction. This expression will be used in the next section 
for writing the manifestly covariant current operators. 

Finally, we recall that, in all the previous eqs.(4.5a-d), the expression of the 
p*Q given in eq.(3.12c) must he used. 

The previous discussion has been focussed on the boost transformation from 
the RF to a GF. However, recalling the general property of eq.(4.3), one 
can immediately verify the equivalence of our model to SF in the case of a 
transformation from a GF to another GF. 



5. The matrix elements of the current operators. Comparison with 
SF. 

In this section we first examine the construction of transition matrix ele- 
ments introducing the RIA; then, we critically discuss the equivalence of our 
formalism with SF. We recall that, in order to compare the theoretical model 
with the experimental data, the electromagnetic and weak form factors can 
be easily extracted from the corresponding current matrix elements [1-3]. 



V'§(P2,P3,V) 



e(P2)e(P3) 
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The main hypothesis of the RIA, as in the nonrelativistic case, consists in 
assuming that, formally, only one constituent quark interacts with the exter- 
nal probe while the others act as spectators. Considering the choice of the 
independent momenta performed in the previous sections, we conveniently 
take the quark #1 as the interacting one and the quarks #2 and #3 as 
spectators. The matrix element calculated according to this hypothesis, is 
then multiplied by a factor 3 to obtain the total amplitude (when considering 
three identical particles). 

In order to construct current transition matrix elements with explicit rela- 
tivistic tensor properties, we shall use the boosted Dirac wave functions of 
eqs.(4.5a-d) and make, in a GF, the integrations over p2 and ps, that are 
the spatial variables of the spectator quarks. According to the impulse ap- 
proximation, these momenta remain unchanged in the initial and final state 
of the scattering process. 

In more detail, we shall denote the four-momentum of the bound system, 
observed in the GF, as Pq. The index G will be set to / and F for the initial 
and final state, respectively. The same notation will be used extensively in 
the following of this section. 

The numerical parameters Vq (introduced in Sect. 3) for boosting the wave 
function from the initial or final RF to the GF, are determined by means of 
eq.(2.3) in the form 

vj^ = P^/Mg (5.1a) 

with 

Mg = ^fP^P^^^ (5.1&) 

As before, the independent components are the spatial ones, i.e. v^. As 
shown in eq.(3.12b), the bound system is in an eigenstate with V = Vg. 
In this work we consider, for dcfinitcncss, elastic transition amplitudes, that 
is with Mg — M, but the method can be generalized to the case of inelastic 

processes. 

For the whole bound system, we introduce the total (measured) four-momentum 
transfer g'^, that is P^ - P^ = g'^ = (g°, q), and = -g^g'* > 0. 

We observe that, on the other hand, the four-momentum (denoted as q'^) 
adquired by the interacting quark ^1, that remains on shell in the scattering 
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process, can be easily calculated from eq.(2.14) and depends on the dynamical 
state of the system. Explicitly, it has the form 

= p'ip - p'i: = < • M/(p2, P3, v^') - < • M/(p2, pa, v,) (5.2) 

At variance with the nonrelativistic impulse approximation, is not equal 
to the measured momentum transfer q^^ [2,3] . 

According to the previous considerations, the current matrix element can be 
written in the following general form 

= 3 j d^p2d^P3d^yd^y' 
V'F (P2, P3, VWf ei f 1 AfilPf{p2, P3, V) 
_L(i + y^y/^SiV - V - Mq)(l + (5.3a) 

= JfMPf -Pi-q) (5.36) 

with 

Jfi = ^ I 4;7T4;rT<^F(P2,P3;vF)[e(p;p)e(p;j,)]^/2 
J e(p2)e(p3) 

Nf ei fi Ni 

[<P;7)e(P;7)]'/V?(P2,P3;V,) (5.3c) 

Let us now comment the previous expressions. 

As anticipated, the factor 3 that multiplies the matrix element, by means of 
the antisymmetry of the wave function, takes into account the contributions 
of the quarks #2 and #3, when these quarks are interacting with the virtual 
photon field. 

The factors in the last line of eq.(5.3a) represent the matrix element, in the 
V representation, of the operator that changes the total momentum of the 

system. 

Eq.(4.5c) has been used to transform eq.(5.3a) into eqs.(5.3b,c). The Dirac 
adjoint wave functions have been introduced multiplying the Hermitic con- 
jugate by 7? O 72 ® 73- 

The factors Hg represent invariant but, in some extent, arbitrary normal- 
ization functions [6] that will be briefly discussed in the following for the 
electromagnetic form factors. 
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The generalized charge operator ci for the interacting quark has been intro- 
duced. The specific form of this operator, in the isospin space, will be given 
in eqs.(5.5) and (5.11) for the electromagnetic and axial current, respectively. 
The symbol Fi denotes the covariant quark interaction vertex. It is given 
by a subset of the 16 Dirac covariant matrices for the quark #1 multiplied 
by spatial functions with definite Lorentz tensor properties. We recall that 
one has the following Dirac matrices: Ti = li, 7^ , 7^7^.... for scalar, vector, 
axial- vector.... matrix elements, respectively. 

The covariant matrix element of the model is Jfj. To clarify the meaning of 
this quantity we recall that, for a single (point-like) spin 1/2 particle, it would 
be represented by the standard expression Jpi — w^^u(P F)eru(P i)wj:j . 

We highlight that our model for the current matrix elements of a compos- 
ite system, represented by eq.(5.3c), is manifestly covariant, according to 
the definition given in the introduction. In fact, in eq.(5.3c) there appear 
covariant integrations over the spectator quark momenta and invariant fac- 
tors. Also, the intrinsic RF wave functions are boosted by means of standard 
Dirac boosts. Finally, the prove of covariance is completed by using standard 
boost transformation properties of the Dirac matrices. In particular: 

S^(v)7°S^(v) =7° (5.4a) 
S^(v)7V5^(v) = 7°L'^^(v)7'^ (5.46) 
and the corresponding transformations for the other Dirac matrices. 

The most relevant case for the study of the hadronic structure is represented 
by the four-vector electromagnetic interaction. In this case one has the fol- 
lowing quark charge operator in the isospin space 

ei = e^^"^ = + ^ (5.5) 

The four- vector vertex can be put in the following phenomenological general 
form 

Ti = = 71^^^ - lar^B • {Kpv^F - Kjv^i) (5.6) 
with the invariant factors 

Fa = Fa{M- P2, P3, , V,) (5.7a) 
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Fb = Fb{M; P2, P3, Vf, Vj) {5.7b) 

Kg = KoiM; vg, P2, P3) (5.7c) 

In the case of a single but nonpoint-like (on shell) particle they represent the 
standard observable form factors 

Fa = F,{Q^), Fb = F^iQ^), Kg = 1 (5.8) 

On the other hand, when considering interacting quarks, Fa, Fb and Kq can 
give an effective representation of all the unknown effects that modify the 
bare quark vertex. Some of this effects can be related to violations of the 
RIA, others to the strong interactions of the constituent quarks. The latter 
are usually interpreted in terms of quark substructure and/or exchange of 
vector mesons between the virtual photon and the quark vertex. 
In any case, at zero momentum tranfer, vertex charge normalization requires 
F4(M; P2, Ps, V(3, Vg) = 1. The second term in eq.(5.6) is related, at zero 
momentum transfer, to the quark anomalous magnetic moment. Further- 
more, in order to represent the dependence of that term on the interacting 
quark momentum transfer, one can take Kq — M/(p2, ps, Vg), obtaining 
as defined in eq.(5.2). 

However, the simplest choice is to consider the interaction of the virtual 
photon with point-like Dirac particles, that is using Fa ~ 1 and Fb = in 
eq.(5.6). In our opinion a relativistic study of the nucleon form factors should 
calculate first these quantities with that choice (by using the nucleon wave 
functions of the quark model), then insert the phenomenological functions 
Fa, Fb and Kq to improve the reproduction of the experimental data. 
The study of the counterterms due to dynamical current conservation and 
the analysis of the contributions due to virtual negative energy states or to 
quark- antiquark pairs can help to construct a more reliable and consistent 
model. 

As for the invariant normalization factors Ag of eqs.(5.3a,c), they can be cho- 
sen considering the requirement of total charge normalization for the matrix 
element at zero momentum transfer. It reads Jq^ = e^^ = +1,0 for the pro- 
ton and the neutron, respectively. This condition is automatically satisfied 
(considering the antisymmetry of the wave function and the normalization 
of the Dirac spinors of eq.(4.1b)) by 
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(5.9) 



Within this theoretical framework, various numerical calculations have been 

performed for the nucleon electromagnetic form factors, by using different 
constituent quark nucleon wave functions. The results, in good agreement 
with the new experimental data, show the essential role of relativity in such 
calculations and the rehabihty of the RIA [2,7] as a starting point for the 
study of the electromagnetic response of the nucleon. 

Similar calculations have been also performed for the study of the axial nu- 
cleon form factor [3,4]. In this case, the quark interaction vertex is taken as 
the axial- vector Dirac matrices 

f 1 = iH (5.10) 
and the axial charge as an isospin raising operator, that is 

d = ei"^ = ri+ (5.11) 

In this case no charge normalization condition can be found and the form 
of the vertex spatial functions and of the normalizations factors should be 
carefully studied [4,6]. 

We now turn to discuss the comparison of our manifestly covariant matrix 
element of eq.(5.3c) with that of the SF [2,6]. For definiteness we refer to 
a four-vector electromagnetic vertex Fi = 7^ with the normalization factors 
given in eq.(5.9). 

To this aim we shall transform our expression of eq.(5.3c) into the SF. We 

divide this procedure into the following three steps. First (i), we obtain the 
rotation matrices of the spectator quarks; second (ii), the rotation matrices of 
the interacting quark; finally (iii), the momentum 5 functions of the spectator 
quarks. 

(i) We now reproduce the rotation matrices of the spectator quarks {i = 2, 3). 
For these quarks, taking into account eqs.(3.8b) and (4.5d), the momenta of 
eq.(3.12c), in eq.(5.3c) one has the following spinorial bilinear quantities 

Sa^'a, = <,«+(p:^)5f (v^)7°i^f (v,)«(p:,)w,, (5.12) 
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By means of eq.(5.4b) one can write 

5f (vF)7°5f (v,) = 7°[5f (v^)]-^Sf (v,) = 7fSf (-v^)Sf (v,) (5.13) 

We now consider the product of the two Dirac boosts in the last equation. 
We recall that the corresponding boosts on the spectator momenta are 

P6(P*/; V/) = Pi (5.14a) 

P6(pIf; vf) = Pi (5.146) 
The last equation can be rewritten as 

P6(Pi; -Vf) = PiF (5.14c) 

In consequence, applying successively (composing) the boosts of eqs.(5.14a) 
and (5.14c), one obtains the following total boost 

P6[P6(P*7; v/) ; -Vf] = p*^ (B.Ud) 

Note that for the Dirac spinors the corresponding boost is the product 
Bf[—\p)Bf{wi) of eq.(5.13). We use for that product, applied to w(p*/) 
the property of eq.(4.3). Then, we insert the result in eq.(5.12). Taking the 
covariant Dirac spinor normalization and the definition of eq.(3.8) for the 
Wigner rotations, one finally obtains 

= <iji?vy(p:,;i^-^(vF)5(v,))] 

= E^;!' [^^(pa>^(vF))]<;jit:,^(p:F,5(v7))] (5.15) 

A. 

The second equality is directly obtained without composing the two succes- 
sive boosts. 



(ii) As for the Wigner rotations of the interacting quark, by means of eq.(4.3) 
and inserting two complete sets of spin states, we introduce the following 
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identity 

Dlil[Rw{plj,B{vj))] (5.16) 

(ill) Let us now consider the spatial integrations over the spectator momenta 
of our eq.(5.3c). We introduce two 5 functions in the following way: 

J £p2d^p3 .... = J d^p',d^p',d^p2d'ps 6ip', - p2)S{p', - Pa) .... (5.17) 

Furthermore, in cq.(5.3), and p'^ are then taken as the arguments of the 
final state wave function. The rest frame final momenta are considered as 
functions of those momenta. The same holds for the initial state, taking p2 
and p3 as arguments. 

We now replace the integration variables p^ (i = 2, 3) and the primed ones 
with the corresponding rest frame momenta p*. For the initial state mo- 
menta, a transformation factor must be introduced according to the following 
equation 

d'p, = 4^ d'p; (5.18) 

An analogous equation holds for the final state primed momenta. We can 
identify 

P* = Pu (5.19a) 

Pf = pIf (5.196) 

and use cqs.(5.14a,b), respectively, to express Pi and p^ as functions of the 
intrinsic momenta. 

Considering eq.(3.9), the rest frame momenta are easily replaced by the Ja- 
cobi momenta as integration variables. 

Collecting all the previous results, our electromagnetic current matrix ele- 
ment is put in the SF, giving 
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<l:'[^H/(p'i*,i?(vi.))]<,«(p;)ei ^^Mvi)wx,D'll^[Rw{vl.B{yi))] 

(^(P2 - P2)5(P3 - Ps) • me(p2)e(p3) 
[e(pr)e(p-)e(p-)e(pt)6(p;)6(p;)]-V2 

V'^S.S(Pp'Pa) (5-20) 
where a sum over the repeated indices is understood. 

We note that the previous expression, that has been shown to be equal to 
eq.(5.3c) for the electromagnetic interaction, is cohincident with eqs.(2),(3) 
and (10) of ref.[6]. Apart from a (probably not relevant) normalization factor, 
our expression is also equivalent to the result of ref . [2] . 

After verifying the equivalence of our covariant matrix element with SF, we 
conclude observing that our expression of eq.(5.3c) presents the following 
advantages with respect to SF. 

(i) As discussed above, it is manifestly covariant. 

(ii) It is more compact, in the sense that it contains only two three-dimensional 
integrations over the spectator momenta with respect to four integrations of 
the SF. 

(iii) Well known Dirac spinors and Dirac boost matrices are used instead of 
rotation matrices of Wigner rotations. 

These features allow for studying the possibility of deriving an expression for 
a dynamically conserved current by means of a suitable procedure of minimal 
coupling substitution. The results of this investigation will be presented in 
subsequent works. 
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